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(N ■ The aim of this paper is to provide an original method of constructing very effective near- 

cloaking structures for the conductivity problem. These new structures are such that their first 
Generalized Polarization Tensors vanish. We show that this in particular significantly enhances 
the cloaking effect. We then present some numerical examples of Generalized Polarization 
I ' Tensors vanishing structures. 
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1 Introduction 

\Q • The central problem in the electrical impedance tomography is to reconstruct the unknown con- 

ductivity distribution of a conductor using the boundary measurements, the Dirichlet-to-Neumann 
on : (DtN) map. Even if unique identifiability by the DtN map holds for wide class of conductivity 

distribution (0 [TS1 E] to cite only a few), Greenleaf et al. [5] found a (singular) conductivity 



distribution whose DtN map is exactly the same as the one associated to the constant conductivity 
distribution. They use a change of variables scheme to create the desired conductivity distribution. 
They push forward the material constant by the transformation blowing up a point to a ball or a 
disk. It turned out that the change of variables scheme can be applied to cloaking: Pendry et al. 
[T6| and Leonhardt [9] used similar idea to initiate the research on cloaking. Cloaking is to make a 
target invisible with respect to probing by electromagnetic waves. Since then, extensive work has 
been produced on cloaking in the context of conductivity and electromagnetism. We refer to [I] 
(also [3]) for recent development on the cloaking. It is worth mentioning that there is yet another 
kind of cloaking in which the cloaking region is outside the cloaking device, for instance, anomalous 
localized resonance [HI [14] . 

The change of variables based cloaking method uses the singular transformation to boost the 
material property so that it makes a cloaking region look like a point to outside measurements. 
However, this transformation induces the singularity of material constants in the transversal di- 
rection (also in the tangential direction in two dimensions), which invokes the difficulty both in 
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the theory and applications. To overcome this weakness, so called 'near cloaking' is naturally 
considered, which is the regularization or the approximation of singular cloaking. In 8 , instead of 
the singular transformation, Kohn et al. use a regular one to push forward the material constant 
in the conductivity equation, in which a small ball (of radius p) is blew up to the cloaking region. 
They estimate that this near-cloaking can be approximated to the perfect one with the order of p d 
in the space of dimension d. 

The purpose of this paper is to propose a new cancelation technique in order to achieve en- 
hanced near-invisibility. Our approach is based on the multi-coating which cancels the generalized 
polarization tensors (GPTs) of the cloaking device. We first design a structure coated around an 
inclusion to have vanishing GPTs of lower orders and show that the order of perturbation due to 
a small inclusion can be reduced significantly. We then obtain near-cloaking structure by pushing 
forward the multi-coated structure around a small object via the usual blow-up transformation. For 
the conductivity equation, we show that the order of near-cloaking is p 2N using N coatings, which 
is a significant improvement over p 2 approximation obtained in [5]. We give numerical examples 
for the material parameters and the thickness of the layers for the GPT vanishing structures. 

This paper is organized as follows. In the next section we derive the multi-polar expansion 
of the solution to the conductivity equation which is slightly different from the usual one, and 
define the (contracted) generalized polarization tensors. In section 3 we characterize the GPT 
vanishing structures. In section 4 we show that the near-cloaking is enhanced (to p 2N ) if the GPT 
vanishing structure is used. In section 5 we present some numerical examples of the GPT vanishing 
structures. We end this paper with a brief conclusion. 

Even though we consider only two dimensional conductivity equation in this paper, the same 
argument can be applied to the equation in three dimensions. The multi-coating technique devel- 
oped in this paper can be applied to the Hclmholtz equation to enhance the near cloaking obtained 
in [HI 13 [TT] . The results for the Helmholtz equation will be presented in the forthcoming paper. 

2 Far-field behavior of the solution 

Let f2 be a domain in R 2 containing possibly with multiple components with Lipschitz boundary. 
For a given harmonic function H in R 2 , consider 

V-(a oX (M 2 \n)+a X (n))Vu = inR 2 , 
u(x) - H(x) = 0(|x| _1 ) as |x| -> oo, 

where Co and a are conductivities (positive constants) of R 2 \ ft and f2, respectively. The solution 
u to (|2.1[) admits the multipolar expansion [1] 

(u-ff)(x) = i -^ r ^r(x)M a ^ J ff(0), |x|->oo, (2.2) 

where M Q( g = M a p{Vi, are the Generalized Polarization Tensors (GPTs) associated with the 
inclusion fi and the conductivity contrast and T(x) is the fundamental solution of the Laplacian, 
i.e., 

r(x) = i-ln|x|. 

Here and throughout this paper a = {a\, 02) and ft = (/3i, P2) are multi-indices and \a\ = ct\ + oti. 

We seek an expression of the multipolar expansion which is slightly different from (|2.2I) . For 
multi-indices a, j3 with a = n, |/3| = n, define (a™) and (frg) by 

^ ««x" = r n cos n9 and ^ b$x? = r n sin n6, (2.3) 

\a\— n |/3|— n 



2 



and define 



M cc = 



M sc = 



E E <%*pM al) , 

\a\—m |/3|— n 

E E <bpM aP , 

\ct\—m |/3|— n 

E E 

\ct\— m \(5\—n 
E E Og^a|9. 



(2.4) 
(2.5) 
(2.6) 
(2.7) 



a\—m \f3\—n 



We call these coefficients the contracted GPTs. 
Using the expansion of In |x — y|, we have 



C— "Hi _i 
E L -j-9xr(x)y" = — 



l«|= 



27rn 



cos n9 



x smnffx . 
r' cos ?w v H — r" sin ?w v 



~,n y 



y y 1 r n ■ y 



n>0, (2. 



where x = r x (cos # x , sin# x ) and y = r y (cos # y , sin# y ), which is valid if |x| 3> 1 and y £ 957, and 
hence 



i — a x r(x) = — 



cos n6 x sin n6> x 
a n = 1- o„ — 



(2.9) 



If the harmonic function H admits the expansion 

oo 

H(x) = H(0) + E r n (a c n {H) cosn0 + a s n (H) sinn<9)) 



n=l 



with x = (r cos 9, r sin 0), then we have 

d a H(0) 



a c n {H)a n a +a s n {H)b n a . 



(2-10) 



Plugging (|2.9p and (|2.10j) into (|2.2p we have the following formula 



cos m# 



(« - h)(x) = e o — -E( M ™>«( ff ) + 



27rmr TI 

m— 1 n— 1 

OO ^ OO 

- E E( M -«<(^) + M ™<(H)) as |x| oo. (2.11) 

m=l n— 1 

We emphasize that (|2.11l) is valid even if er (the conductivity of f2) is not a constant but a 
variable. In the next section we will design a conductivity distribution a of f2 so that M^f„ = 
M^f n = M^ n = M^ n = for all m, n < N for a given integer N. We call such a conductivity 
distribution GPT -vanishing structure or coating of order iV. 



3 GPT vanishing structures 

To obtain GPT-vanishing structures, we use a disc with multiple coatings. The idea comes from 
Hashin's neutral inclusion which is a disc with a single coating [12 . The special property of 
the neutral inclusion is that it does not perturb the uniform fields outside the inclusion, which 
is equivalent to the first order polarization tensors of the inclusion vanishing. In other words, 
Hashine's neutral inclusion is a GPT-vanishing structure of order 1. 
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Let n be a disk of radius r±. For a positive integer TV, let < rjv+i < < ■ ■ ■ < ri and define 

A, := fa+i < r < r,-}, j = 1,2,..., iV. (3.1) 

Let Ao = R 2 \f2 and A^ + i = {r < rpj + i}. Set crj to be the conductivity of Aj for j = 1,2,..., N+1, 
and (To = 1- Let 

AT+1 

Let M^jjcr], etc, denote the (contracted) GPTs associated with a (and ft). Because of the sym- 
metry of the disc, one can easily see that 

AC>] =M s ^ n [a] =0 for all m,n, (3.3) 

M% n [o) = AC>] = ifm^n, (3.4) 

and 

M-[<7]=M-[<7] for all n. (3.5) 

Let M n = M^ n , n = 1, 2, . . ., for the simplicity of notation. 
To compute , we look for solutions Uk to 

V • dVu = in M 2 (3.6) 

of the form 

b (k) 

u fe (x) = a^V fc cosfc6> + -jj-coskd in Aj, j = 0, 1, . . . , N + 1, (3.7) 
with a fe ^ = 1 and = 0. Then is the solution to (|2.ip with if(x) = r fc cosk8, and satisfies 

(u-#)(x) = -^-cosk9 as|x|^oo. (3.8) 

Hence, we have 

M k = -27r/c6 fc) . (3.9) 
We observe that the solution to (|3.6[) is given by 



Thus 



Ufc (re ie )= (r*-— ^jcosfcfl in A,. (3.10) 



since otherwise Ufc(x) = or Ufe(x) = constant in |x| < r for some r > 2. Thus we have 



rfe Mfc > Q and rk + M fc > q for any r > 2 , 

2i:kr k 2irkr k 



and hence 

|M fe | < 27rfc2 2fc for all fc. (3.11) 
The transmission conditions on the interface {r = rj} are given by 



7 ( fe ) r fe 













i 



















(3.12) 



*J 4^ - I = <*-i I "SVT " « I ■ (3-13) 
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Thus we have 



and hence 



" (fc)~ 


1 




" 2^ . 



U j +(7j-l 

(aj - a, ,)/•;'' 



'3 

(J j +CTJ-! 



„-2fc 





" (fe) " 
oJJ x 




6 (fe) 



(fe) 

JV+1 





JV+l 



.2A- 



-2 ft 





' 1 " 




b (k) 



Let 



p (fe) 



- (fe)' 

Pll Pl2 

(fe) (fe) 

P21 P22 



N+l 



Then, we have 



(J.; + tr 7 _i 



2/,- 



(Tj + (Tj-l 



(fe) 
P21 

(fe) ' 
P22 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



Note that Mk = if and only if p( fc ) is an upper triangular matrix, i.e., = 0. Let 

,N. (3.18) 



A, .— , J — 1, 



We now characterize GPT-vanishing structures. 

Proposition 3.1. If there are non-zero constants X±, . . . , Aat + i (\Xj\ < 1) and r% > . . . > r/v+i > 
such that 



N+l 

n 



1 

NT 3 



2k 



is an upper triangular matrix for k = 1, 2, . . . , N, 



(3.19) 



£/ien (f2, <r), given by H3.1\) . HS.2\) . and 113. 18)) . is a GPT-vanishing structure with Mk = for k < N . 
More generally, if that there are non-zero constants X%, A2, A3, . . . (\Xj\ < I) and r\ > ri > r% > . . . 
such that v n convevges to <x positive Huuibev, say > 0, and 



n 

3=1 



1 



a, 



„2A- 



2/v 



is an upper triangular matrix for every fc, 



(3.20) 



i/ien (fi, cr) 7 given by \3.1\) . i3.2\) , and \3.18)) . is a GPT-vanishing structure with Mk = for all k. 

Note that (|3.19p and (|3.20l) are nonlinear equations. For example, if N = 3, equation (|3.19l) is 
simple and reduces to 



Xirf + A 2 rf + A 3 rf + X l X 2 X 3 ri k r^rf = 0, k = 1, 2, 3 



(3.21) 



It is quite easy to show that it admits infinitely many solutions Ai, A2, A3 (|Aj| < 1) and r± > r-i > 
r3 > 0. However, as N gets larger, solving analytically equation (|3 . 19[) seems too complicated, 
and even proving existence of solutions to (|3.19l) or (|3.20p seems to be quite challenging. We 
present a simple numerical method to find the GPT-vanishing structures in Section [5J which is 
also important from a practical point of view. These numerical evidences show us that (|3.19l) has 
solutions, even though we are not able to prove it. 
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4 Near-cloaking using GPT-vanishing structures 



In this section we achieve enhanced near-cloaking by using GPT-vanishing structures. 

Let B r be the disk centered at the origin with radius r and B = B\. Let Aj, j = 0, 1, . . . , N+ 1 
be defined by (|3.1[) with r\ — 2 and rjy+i > 1- Let a be the conductivity distribution defined by 
p.2p . For a given domain f2, we denote the DtN map of f2 with the conductivity a as An[c], which 
is given by 



AoM(/) = f 7— 



where it is the solution to 



V ■ crVit = 



in 17, 
on 951. 



(4.1) 



(4.2) 



Proposition 4.1. Let a be a conductivity profile on R 2 defined as VS. 2^1 and for a small positive 
constant let 



*i(x) = -x, xeK 2 . 
P 



(4.3) 



2ks- l p 2k M k [a] e±lfce _ 



Then the following holds for k = 0, 1, 2, . . . arte? s > 

(A S .[-^]-A B8 [l])(e^) = ^_ Mfe[(r]p2fc 

Before proving Proposition ^. 1[ let us make a few remarks. If cr is a GPT-vanishing structure 
of order N, i.e., M k = for all k < N, then (14.41) shows that 



(4.4) 



Ab b [a o * i ] - A Bs [1] {e ±lk& ) = for all |fc| < N 



(4.5) 



In other words, As[cro vpi] and As[l] cannot be distinguished by slowly oscillating Dirichlet data. 
Moreover, the complete GPT-vanishing structure (|3.20p is achieved, then 



A Bs [cro*i] - A Bs [l])(e 



for all k, 



(4.6) 



which would yield the perfect cloaking. 

Proof of Proposition \4-l\ If u is the solution to 

V • (ffofi)Vu = 

P 

u = f 



in B. 



on |x| = s, 



then u := uo ^ p satisfies 



Moreover, we have 



V • crVu = in B± , 

P 

v w = /o* p on |y| = 



(4.7) 



(4.8) 



du 
dv 



„ y) = P^r 
|y|=f 01/ 



|x|= S 



(py)- 



Therefore, 



A B , M(/ o 4< p ) = /9 A Bs [cr o * i ](/) o * p . 

P P 

To compute [a], set fie 16 ) = e ife<9 for a fixed fc € N (fc > 0). Then, the solution u k to 

p 

6 (fe) 

fi fc (x) = afr k j ke + -±-e ike in A h j = 0, 1, . . . , N + 1, 



(4.9) 



is given by 



G 



with 



a^s k p- k + bi k) s- k p k = l 



(4.10) 



and bP +1 = 0. From ([3J4]l and (|3X6j> , we get 



,(*0 




(fc) (fc) 

Ml Pl2 

(fc) (fc) 

P21 P22 



7 W 

L(fe) 



and hence 



,(fe) _ P21 (k) _ M fc [tr] w 

B W a ° " 27rfc ■ 
P22 



Substituting it into (|4.10p . we have 



-M k [a}p k s k 



b (k) = - 

2nks 2k - M k [a}p 2k ' 



and we obtain 



9j/ 



|y|=1 



jt£(fli*),y*_^),-y) e iM 



fe^(4 fc) S V- fe + ^ fe) s-V)e ife9 
s u s fc+1 



^ fc0 o sfe+ i e 



From (|4.9[) we have 



Ab.[1]( 



ifce \ 



2fe-V 2fc M fe [ ( r] 



□ 



2vrfcs 2fe - M k \a\p 2k 

Thus we get (|4.4p . The same argument works for fc < as well. The proof is complete. 
If / admits Fourier expansion /(e ) = X^fcez /fc e > then we have from (14.41) 

^°. i i-A.[ii)(/ )= g jffr^t 1 H /t '" 

For the GPT-vanishing structure of order TV, we have M k = for all fc < N and from (13. lip that 

p 2k \M k \ < 2irk(2p) 2k for all k > N. 

Since p <C 1, we obtain 



p 2k \M k \ < Cp 2N+2 forallfceN 
for some constant C independent of fc, and hence 

A B >o¥ i ]-A B .[l] <C P 2N+2 , 



(4.12) 



(4.13) 



where the norm is the operator norm of H^ipBs) into #-V2(aB a ). Here H 1 / 2 {dB s ) is the usual 
Sobolev space of order 1/2 on dB s and H~ 1 / 2 {dB s ) is its dual. 
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We mention that the solutions for GPT- vanishing structure exist numerically even when er/v+i — 
0, which is equivalent to prescribe the insulating condition on the boundary of the inner core. Some 
examples with ctat+i = and rj\r+i = 1 are given in the following section. In such structures, the 

conductivity a o $i is in |x| < p, fluctuates in p < |x| < 2p, and is 1 in |x| > 2p. In order to 

p 

have near-cloaking device in 1 < |x| < 2 with the core part |x| < 1 insulated as was considered in 
[5] , one may use transformation as was done in [H [TH] . 

For a given domain Q and a subdomain B C f2, we denote the DtN map of fi with the 
conductivity a as Af2,s[cr], which is given by 



An,sM(/) = o 



du 
dv 



an 



where u is the solution to 



( V • aVu = 

8u =Q 
ov 



in n \ B, 
on dB, 
on Oft. 



By (|4.4I) . we have 



(a B2 , b >o^]-a B2 , [i])(/) = ]T 



|fc|(f) 2|fc| M| 



h-kZ 



27r|fe| 



and hence by (|4.12l) 

AB !A h*i] - A fla ,0[l] 
We define the transformation F p : B2 — > B2 by 



< Cp 



. 2|fej 



2W+2 



Mi, 



-/fee 



-3-3p 



1 



2p 3-2p 



for § < |x| < 2, 
for p < |x| < §, 

for |x| < p. 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



Then one can easily see that 

Ai3 2 , Sp [CTO*i] = A Bq , Bi [(F p )*(cro1<i) 
p L p 

where 



H)H*i)(w P ) J 



\det(DF p )\ 

We then obtain the following theorem, which is the main result of this paper, from (|4.16|> . 
Theorem 4.2. There exists a constant C independent of p such that 

Lh j: L (F,Maoti)| -As 2)0 [l] <C P 2N + 2 . 



A 



(4.18) 



(4.19) 



5 Numerical examples 

In this section we present some numerical examples of ci, . . . , ctat+i and r\ > . . . > fAr + i satisfying 

dun. 
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We fix N and Tj = 2 — ■2— L for j = 1, . . . , N + 1. We then solve the following equation for 

<T = (<Tl, . . . , OjV+l) 

<r >->• M fe [<r] = for fc= 1,...,JV. (5.1) 

Since (|5.ip is a nonlinear equation, we solve it iteratively. Initially, is set to be 2^ _1 ^, 



j = 1, . . . , N + 1. We iteratively modify er^' = (er^ , . . . , as 

a (i+i) =o .W_At b (0 ) 



r (0 



where is the pseudoinverse of 



A; := 



cKM!,...^) 



and 



b« = 



dc 



Mi 



Example 1. Figure I5TT1 shows computational results of the conductivity a for N — 3,6,9. It 
clearly shows that the larger N is, the more a fluctuates. One interesting thing to observe is that 
ctjv+i takes values of 1.9695, 0.9791, 1.0029 for N — 3, 6, 9, respectively: they are getting closer to 
1 which is the conductivity of the exterior part. 

Example 2. Here we see what happens if the conductivity ctjv+i of the core D is fixed with a value 
different from that of exterior, namely 1. We set N (the number of layers) to be 9 and trio = 5 
and 0.2. The numerical results are illustrated in Figure [5~^1 Conductivity profiles fluctuate more 
drastically and takes greater values than those in Example 1. 

Example 3. This example is for the near-cloaking for which the boundary of the core (r = 1) is 
insulated, and hence the conductivity of the core is set to be 0. Figure [5731 shows the results of 
computation when TV = 3, 6: the conductivity fluctuates on coatings near the core. When N = 3, 
the maximal conductivity is 5.5158 and the minimum conductivity is 0.4264. When N = 6, they 
are 11.6836 and 0.1706. 



6 Conclusion 

We have obtained new near-cloaking examples for the conductivity problem. We have shown that 
the GPT-vanishing structures can be used to enhance the near-cloaking. The GPTs up to the 
N-th order can be canceled using (N + 1) layers with different conductivity parameters. To make 
the numerical procedure simple, we have assumed that the layers are concentric disks centered at 
the origin with specific radii. The numerical simulations show that the GPT-vanishing structure 
exists with these radii. As mentioned in Introduction, the multi-coating technique can be applied 
to enhance near-cloaking for the Helmholtz equation. 
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